Stability of complex networks under the evolution of attack and repair 
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With a simple attack and repair evolution model, we investigate and compare the stability of 
the Erdos-Renyi random graphs (RG) and Barabasi- Albert scale-free (SF) networks. We introduce 
a new quantity, invulnerability I(s), to describe the stability of the system. We find that both 
RG and SF networks can evolve to a stationary state. The stationary value 7 C has a power-law 
dependence on the average degree {k) rg for RG networks; and an exponential relationship with 
the repair probability p s f for SF networks. We also discuss the topological changes of RG and SF 
networks between the initial and stationary states. We observe that the networks in the stationary 
state have smaller average degree (fc) but larger clustering coefficient C and stronger assortativity 
r. 

PACS numbers: 89.75.Hc, 87.23.Kg, 89.75.Fb 



I. INTRODUCTION 

Complex networks are an essential part of modern soci- 
ety. Many social 0], biological @, transportation [J0| 
and communication i systems can be cast into the form 
of complex networks, a set of nodes joined together by 
links indicating interactions. The investigation of com- 
plex networks was initiated by Erdos and Renyi pj in 
1950s. They introduced a kind of simplest and most 
straightforward complex networks, described as random 
graphs. In a random graph, N nodes are connected ran- 
domly with the probability p rg . This model has guided 
our thinking about complex networks for decades since 
its introduction. In 1999 Barabasi and Albert 0, 0, 0] 
began to put the emphasis on the network dynamics and 
constructed networks with the algorithm of growth and 
preferential attachment. According to the BA model, the 
network grows over time by the addition of new nodes 
and links. A node newly added to the network randomly 
selects m other nodes to establish new links, with a selec- 
tion probability that increases with the number of links 
of the selected node. One of the most relevant is given by 
the scale-free nature of the degree distribution P(k), de- 
fined as the probability that a randomly chosen node has 
degree k (i.e., it is connected to other k nodes). In math- 
ematical terms, the scale-free property translates into a 
power-law function of the form P(k) ~ k 1 . 

Recently the security of complex networks to the ran- 
dom failures or intentional attacks has attracted a great 
deal of attention [Til Il2l Il3l| . The random failure is sim- 
ulated as the deletion of network nodes or links chosen 
at random, while intentional attack as the targeted re- 
moval of a specific class of nodes or links. The work by 
Albert et al demonstrated that scale-free networks, 
at variance of random graphs, are robust against random 
failures of nodes but vulnerable to intentional attacks. 
The analytical approach 0, 0] has also been developed 
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for investigating the error and attack tolerance. 

However, despite the fine work of studies on network 
tolerance, little effort has been made on the repair of the 
attacked networks. In this paper, we propose a simple 
model and study the evolution of complex networks un- 
der continuous attacks and repairs. This effect can be 
pictured for the organism: the organism has to expe- 
rience various DNA damages and corresponding repairs 
during its growth. The effect can be also exemplified by 
the Internet: the development of Internet always follows 
the contest of virus attacks and protections. Our model is 
aiming at investigating the stability appearing in the evo- 
lution of complex systems frequently attacked by internal 
or environmental agents. The model shows that complex 
networks will become more stable in the evolution and 
eventually reach the stationary states. We also discuss 
the topological changes of complex networks before and 
after the evolution. It will be useful for understanding 
the dynamics of complex systems. 

II. EVOLUTION MODEL 

We start by constructing a network according to the 
Erdos-Renyi random graph (RG) model or Barabasi- 
Albert scale-free (SF) model. The algorithms of RG and 
SF networks have been mentioned in the Introduction. 
The dynamics of our model is defined in terms of the 
following two operations: 

1. Attack: Find a node with the maximum degree 
kmax and remove all its links. (If several nodes 
happen to have the same highest degree of connec- 
tion, we randomly choose one of them.) 

2. Repair: Reconnect this node with the other nodes 
in the network with probability p re . 

Then, the evolution comes into the next Monte Carlo 
time step. The random repair is considered because most 
of the time our operation to the system is more or less 
blinded due to incomplete information. In addition, the 
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FIG. 1: Plots of kmax versus time step s with size N = 10 
000 and p re = 0.005 for (a) random graphs; (b) scale-free 
networks. 



FIG. 2: Plots of I(s) versus step s(million) with N = 10 000 
and pre = 0.005. 



size of the system remains unchanged during the evolu- 
tion since the way of simulated attacks is the deletion of 
links. In order to reduce the number of parameters in 
the model, we set the repair probability p re the same as 
the connection probability p rg in RG networks. This set 
also implies that the information does not increase from 
initial construction to later repairing. 



III. RESULTS AND DISCUSSIONS 

A. Stability of complex networks 

We have performed extensive numerical simulations to 
study the evolution of RG and SF networks. It can be 
easily seen that the maximum degree k max has a ten- 
dency to decrease in the evolution. In Fig. 1, we give a 
snapshot of the maximum degree k max versus time step s 
with N — 10 000 nodes and repair probability p re — 0.005 
for RG and SF networks. We find that the fluctuation 
of k max is very large for RG networks. For SF networks, 
kmax decreases very steeply at first and then steadily 
with s. The behavior of scale- free networks is rooted in 
their extremely inhomogeneity in which the removal of 
a few highly connected nodes dramatically reduces the 
connections of the networks. 

From intuition, a node with less links to others will 
be attacked less frequently. Thus a network with smaller 
maximum connection degree is more stable. To describe 
the stability of the network, we introduce a new quality, 
Invulnerability I(s). The definition of I(s) is analogous 
to that of the Gap G(s) in the Bak-Sneppen evolution 
model [l7L ITsI ] , which is used to simulate the evolution of 
species. Considering an evolution of network with maxi- 
mum degrees k max (l), k max (2),..., k max (s), invulnerabil- 



ity I(s) at time s is defined 

I{s) = l/Miri{U(t)} for i < s , (1) 

i.e., the inverse of the minimum of all the maximum de- 
grees in the evolution before moment s. Initial value of 
I(s) is equal to l/k max (l). I(s) reflects the attack tol- 
erance of the network. When I(s) is small, the network 
is vulnerable to attack and behaves less stability. Ob- 
viously from definition, I(s) is non-decreasing function 
of evolution, and some fluctuations in k max have been 
filtered out. 

In Fig. 2 we plot the evolution J(s) versus s(million) of 
RG and SF networks under N = 10 000 and p re = 0.005. 
We observe that in both networks J(s) increases very 
quickly at small s but slowly at large s and finally reaches 
a stationary value I c when s is large enough. The in- 
crease of I(s) indicates that the system is getting more 
and more stable in general by experiencing continuous 
attacks and repairs. This figure is similar to the enve- 
lope function of Bak-Sneppen evolution model. Without 
interference from outside world, the network evolves to a 
stationary state, and the process takes place over a very 
long transient period. 



B. Stationary states of random graphs 

In this section we will present the properties of sta- 
tionary value V c 9 in the stationary state for RG net- 
works. In RG networks, the average degree (k) rg is deter- 
mined by the connection probability p rg and the size N, 
(k) r g = PrgN. With (k) rg = 20 fixed and p re = p rg , we 
find that I^ 3 stays unchanged at 0.2 with size N ranging 
from 100 to 10 000. This result is interesting, because it 
shows that the stationary value I^ 9 depends not on the 
network size N, but the average degree (k) rg . 
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FIG. 3: Plot of Ic 3 versus {k) rg with size N ranging from 100 
to 10 000. 
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FIG. 4: Plot of I*f versus repair probability p a f with size 
N = 1000 (circles); 5000 (squares); and 10 000 (triangles). 



To get the relationship between l r c 9 and (k) rg , we show 
in Fig. 3 I^ 9 as a function of (k) rg in a log-log plot. We 
find that l r c 9 has a power-law dependence on the average 
degree {k) rg , 



i r c s ((k) rg ) ex <fc)-;, (2) 

where the exponent r is about 1.5. Fig. 3 illustrates that 
after the network has relaxed to the stationary state, the 
stability of the network will increase rapidly with the 
decrease of average degree (k) rg . Thus, when the (k) rg 
is small, i.e., less communications and interactions in the 
network, the system will be more stable. 



C. Stationary states of scale- free networks 

The properties of stationary value I** in the stationary 
state for SF networks will be illustrated in this section. 
In Fig. 4, we plot J** as a function of repair probability 
p s f in a semi- log coordinate with N — 1000 (circles); 
5000 (squares) and 10 000 (triangles). We find that I S J 
decreases with the increase of N or p s f. While large N or 
p s f means the large number of links in networks. In other 
words, the stability of the system relies on the internal 
interactions within the system, that is, the stronger the 
interactions, the less stable the system. From Fig. 4 we 
can also observe that in this semi-log plot 1^ has an 
approximately linear decrease with p s f under fixed N, 
i.e., an exponential relationship between I** and p s f, 



I s J(p sf ) ~ eM-a(N) Psf + b(N)), (3) 

where a(N) and b(N) are two fitting parameters related 
to the size N. 



TABLE I: Comparison of topological changes in initial and 
stationary states between RG and SF networks under N — 10 
000 and p re = 0.005. 
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D. Topological changes of complex networks 

In order to offer further information of complex net- 
works in the stationary states, we compare some struc- 
tural properties of the networks in the initial states with 
those in the stationary states. In Table [IJ we display the 
topological changes of RG and SF networks under N — 
10 000 and p re = 0.005. We use the superscript i and s 
to distinguish the initial and stationary states. 

Average degree. We find that the average degree in 
stationary states (k) s decreases dramatically for RG net- 
works but slightly for SF networks. This result is rooted 
in the homogeneity of RG networks. Since in RG net- 
works all nodes have approximately the same number of 
links, majority of nodes in the network will be impacted 
by the attacks and, hence, the rapid decrease of (k). 

Characteristic path length. The characteristic path 
length L is the average length of the shortest paths be- 
tween any two nodes in the network. L s in the stationary 
state is a little larger than L l , which exhibits that it will 
take a little longer for the nodes in the network to com- 
municate with each other after the evolution. 

Clustering coefficient. The clustering coefficient of 
node i is defined as the existing numbers rij among 
the links of node i over all the possible links, that is, 
Ci = rii/[ki(ki — l)/2]. The clustering coefficient of the 
whole network is C = C s for RG and SF net- 

works are increased by 99% and 72%, compared with C l 
in the initial states. It suggests that the nodes are highly 
clustered in the stationary states. 

Assortativity. The assortativity r in the range — 1 < 
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FIG. 5: A comparison of degree distributions under N = 10 
000 and p re = 0.005 in initial (circles) and stationary (aster- 
isks) states for (a) random graphs; and (b) scale-free networks 
. Solid lines are the fit curves. 



r < 1 is another interesting feature of complex networks. 
A network is said to show assortative, r > 0, if the high- 
degree nodes in the network preferentially connect to 
other high-degree nodes. By contrast, a network is disas- 
sortative, r < 0, if the high-degree nodes tend to connect 
to other low-degree nodes. The values of r % and r 8 have 
shown that both RG and SF networks turn to strongly 
assortative networks in the stationary states from initial 
weakly assortative or even disassortative ones. 

Degree distribution. Fig. 4(a) and 4(b) display the 
degree distributions P(k) of RG and SF networks in 
the initial and stationary states under N = 10 000 and 
p re = 0.005. For RG networks, the degree distribution 
P{k) s rg in the stationary state is a Gaussian one with 
mean value (k) rg ~ 19 and standard deviation a s rg ~ 7, 
while the P(k) rg in the initial state is a Poisson distri- 
bution with (k) rg ~ 50 and a Tg = ^ (k) rg ~ y/50. For 
SF networks, P{k)\j and P{k) s s ^ are both power-laws, 
P(k) s f <~ k~f , with the similar exponent about 2.4. The 
differences between P(k)lj and -P(fc)s/ are that the fat 
tail in P{k) l a ^ vanishes; and the maximum degree de- 
creases from the initial value of 241 to 50 in stationary 
state. Besides, Fig. 4(a) and 4(b) show that the shape of 
the degree distributions is almost the same for RG and 
SF networks in initial and stationary states. It suggests 
that the general structures of the systems do not change 
even after a long time evolution to the stationary states. 



In general, the system in the stationary state has 
smaller average degree (fc) but larger clustering coeffi- 
cient C and stronger assortativity r. It implies that the 
system after undergoing a long period of evolution is com- 
posed of small highly clustered clusterings, which well fits 
the natural world. In ecosystems, agents autonomously 
form communities, in which individuals arc highly cor- 
related. Commercial airlines, for reasons of efficiency, 
prefer to have a small number of hubs where all routes 
would connect. 



IV. CONCLUSIONS 

In this paper we study the stability of random graphs 
and scale-free networks with a simple attack and repair 
evolution model. We introduce a new quality, invulner- 
ability I(s), to describe the stability of the networks. A 
stationary state with fixed I(s) is observed during the 
evolution of the network. The stationary value of invul- 
nerability I c is found to have a power-law dependence 
on (k) rg for RG networks; and an exponential relation- 
ship with repair probability p s f for SF networks. We 
give further information on the evolution of the networks 
by comparing the topological changes between the initial 
and stationary states. We conclude that the networks in 
the stationary states are composed of small highly clus- 
tered clusterings, which well fits the evolution of complex 
systems in natural world. 

Still, there are many issues to be addressed, such as the 
correlations and the fluctuations in k max in the evolution, 
especially after the stationary state. These fluctuations 
may tell us more about the nature of the stationary state. 
The behavior of the average degree in the evolution also 
may shed some light on the dynamics of the networks. In 
addition, it is worthwhile to investigate whether the evo- 
lution of networks demonstrate self-organized criticality 
(SOC), according to the similarity between the evolution 
of invulnerability I(s) with that of the envelope function 
G(s) in Bak-Sneppen model. All these topics can not be 
covered in this paper, and will be discussed later. 
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